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Abstract 

This paper describes the numerical solution of the problem in classical celestial navigation: 
calculate the position by two lines of position. 

Two cases are treated: the Fix, where the two observations may be considered as 
simultaneous, and the Running Fix, where in general due to the time elapsed between the two 
measurements with the sextant, you must move the first line to the instant of the second one. 


© Andres Ruiz, 1999 
Navigational Algorithms 
San Sebastian - Donostia 
43° 19’N 002°W 
Current version: 2007 
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Variables _ 

P Intercept of a sight 

p = Ho - He 
Towards = + / Away = - 
Z Azimuth (true) 

Measured clockwise around the horizon from 
0° to 360°, is the arc of the horizon between 
the meridian of a place and the vertical circle 
passing through a celestial body. 

B Latitude 

N (+) / S (-) 

L Longitude 

E (+) / W (-) 

Equation of the Line of Position 

The equation of a line in the Cartesian 
plane is given by: 

y = f(x) = m x + b 

Where m is the slope and b is the 
ordinate at the origin. 



From the figure it follows the relationship 
of m y b, with the parameters characterizing 
the line of position, LoP: p and Z 

m = tan(180 9 -Z) = -tan Z 
b = p/cos Z 


Substituting these values into the 
equation and rearranging terms, leads to the 
equation of the LoP: 

p = x sin Z + y cos Z 

Fix by two LoP 

Position (B, L) can be calculated 
analytically as the intersection of two LoP, 
taking the estimated position as the origin of 
Cartesian coordinates. The process merely 
mathematically formulate the graphic 
method used in the nautical chart. 

If both are obtained simultaneously: 

pi = x sin Z1 + y cos Z1 
p2 = x sin Z2 + y cos Z2 


Solving this system of two equations the 
point (x, y) is obtained; the cut between the 
two LoP. The latitude and longitude will be 
obtained from the estimated position (Be, 
Le). 

If p is in nautical miles: 

B = Be + y/60 
L = Le + x/60/COS(B) 

The calculation process is detailed in the 
appendix. The required data is: 

• Estimated Posicion: Be, Le 

• LoPI: pi, Z1 

• LoP2: p2, Z2 

Running Fix 

In the more general case the two LoPs 
will be obtained for different moments in 
time: tl and t2. 

Let R and d be the course and distance 
sailed between tl and t2. The first LoP, 
RA1 (tl), must be moved to the instant of the 
second observation; RA1(t2). In this way the 
intersection of RA1(t2) and RA2(t2) give us 
the position. 

The distance sailed can be calculated as 
a function of time elapsed between the 
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observations and the speed sailed 
constant course: 


d = V (t2-t1) 


at 



Move of a LoP. 

In the figure, K is the distance to be 
moving the first LoP. 


cos a = K/d 
Z1 + a = R 
K = d cos(R-ZI) 

So the intercept was corrected by: 
P 1(t2) = pl(tl) + K 

p1(t2) = pl(tl) + d cos(R-ZI) 


Of course, if both observations are 
simultaneous, tl = t2, the distance is d = 0, 
and the intercept does not change. 


Navigational Algorithms 


A1. Algorithm 

Situacion por dos rectas de altura 
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/ Observacion 1 

UTIi Hs 1p Z-i 

Observacion 2 

UT1 2 Hs 2i Z 2 


A 



Determinant© RA 




r m 



( £, () 


Situacion Observada 

B = Be + y/60 
L = Le + x/6G/COS{B) 


RA 1 


RA 2 

pl.ZI 


p2, Z2 
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A2. Example 
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A3. Software 



Available at the Navigational Algorithms web page. 
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A4. Source code 


/* 

FILE: 2RA.C 

Posicion por 2 Rectas de Altura 

This file contains proprietary information of Andres Ruiz Gonzalez 
Andres Ruiz. San Sebastian - Donostia. Gipuzkoa 
Copyright (c) 1998 
*/ 

#include <math.h> 

#include "../Nav/angulos.hpp" 


void Posicion2LOPSimultaneas ( double Lat 

double *B, 


{ 

double x, y; // [nm] 


double 

double 


Lon, 
L ) 


double pi, 


double Zl, 


double p2, 


double Z2, 


*B = *L = 0; 

double a = SIN(Zl); 
double b = COS(Zl); 
double c = SIN(Z2); 
double d = COS(Z2); 

x = (pl*d-p2*b)/(a*d-b*c); 
y= (p2*a-pl*c)/(a*d-b*c); 

*B = Lat + y/60.0; 

*L = Lon + x/60.O/COS(*B); 

if ( *L > +180.0 ) *L -= 360.0; 
if ( *L < -180.0 ) *L += 360.0; 


void Posicion2LOPNoSimultaneas ( double Lat, double Lon, double pi, double Zl, double p2, double Z2, 

double R, double D, 
double *B, double *L ) 

{ 

Posicion2LOP( Lat, Lon, pl+D*COS(R-Zl), Zl, p2, Z2, B, L ); 

} 


A5. References 

• http://es.wikipedia.org/wiki/Recta 
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